Observation of standard spin-switch effects in F/S/F trilayers with a strong 
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We have measured the superconducting transition temperature T c of F/S/F trilayers using Permal- 
loy (Py=Nig4Fei6) as a strongly polarized ferromagnetic material. For a parallel (P) or anti-parallel 
(AP) alignment of the magnetization directions of the outer ferromagnets, we observe a T c difference 
as large as 20 mK, with a stronger suppression of superconductivity in the P state than in the AP 
state. This behavior is opposite to the recent observations of Rusanov et al., Phys. Rev. B 73, 
060505 (2006) in Py/Nb/Py trilayers, but is consistent with earlier results on trilayers with Ni or 
CuNi alloy as the ferromagnetic material. 

PACS numbers: 74.45.+C, 85.75.-d, 85.25.-j, 73.43.Qt 



The presence of a ferromagnetic (F) material in con- 
tact with a conventional superconductor (S) results in a 
strong mutual influenced The superconducting correla- 
tions penetrate into the ferromagnet and oscillate in sign 
over a very short distance, due to the large energy dif- 
ference between the majority and minority spin bands in 
the ferromagnet. Bilayers, trilayers, and multilayers of S 
and F materials exhibit a wide variety of novel phenom- 
ena, including oscillations of the superconducting critical 
temperatur o 2 i?i 4 and density of states^, and Josephson 
junctions with a 7r-shifted ground state. 6 

In this paper we focus on the so-called "superconduct- 
ing spin switch" first discussed in 1966 by deGennes 7 
and rediscovered in 1999 by Tagiro\£ and by Buzdin, 
Vedyayev and Ryzhanova. 9 Those authors predicted that 
the critical temperature, T c , of a F/S/F trilayer should 
depend on the relative magnetization direction of the 
two F layers, with the smallest T c occurring in the par- 
allel (P) state and the largest T c in the antiparallel 
(AP) state. Those predictions were verified long ago 
by Deutscher and Meunier^i " and more recently by Gu 
et aZ.jii Potenza and Marrows^ and Moraru et ali& 
in a variety of F/S/F systems. It came as a surprise, 
therefore, when Rusanov et alM* recently reported ob- 
servation of the inverse spin switch effect in a series of 
Py/Nb/Py trilayer samples. Although the difference in 
T c between the P and AP magnetization configurations 
was small in that work, the data showed clearly that the 
resistance in the transition region was higher for the AP 
configuration than for the P one. In fact, similar be- 
havior had previously been observed by Pena et aliU in 
F/S/F trilayers made from superconducting YBa2Cu3C>4 
and ferromagnetic Lao.7Cao.3Mn03, with a spin polar- 
ization expected to be close to 100%. Those authors 
interpreted their observations as arising from enhanced 
reflection of spin-polarized quasiparticles at the F/S in- 
terfaces in the AP state leading to a stronger suppression 
of superconductivity, 16 and Rusanov et ali^ claimed that 
the inverse spin-switch behavior is generic for F/S/F tri- 
layers with strong ferromagnets. We believe that the 
mechanism based on reflection of quasiparticles at the 
S/F interface 1 ^ can explain changes in resistance under 



nonequilibrium conditions, but cannot explain differences 
in the equilibrium T c between the P and AP states. Given 
our earlier work showing standard spin-switch behavior 
in Ni/Nb/Ni trilayers^ we were motivated to carry out 
independent measurements of T c in Py/Nb/Py trilayers. 

A series of Py(8)/Nb(d s )/Py(8)/Fe 50 Mn 5 o(8)/Nb(2) 
multilayers (all thicknesses are in nm) was fabricated, 
with thicknesses for the superconducting layer, d s , vary- 
ing between 20 and 150 nm. The samples were grown 
directly onto Si substrates by magnetically-enhanced tri- 
ode dc sputtering in a high vacuum chamber with a base 
pressure in the low 10~ 8 Torr and an Ar pressure of 
2.0T0~ 3 Torr. The thickness of the Py layers were chosen 
to be much longer than the dirty-limit coherence length, 

\J HDp I E ex , which we estimate to be 1.2 nm for the ma- 
jority spin band using E ex — 0.135 eVji£ D F = v F V F /3, 
v F = 0.22 ■ 10 6 m/s and l F = 4 iim. 1? The FeMn layer 
fixes the direction of the top Py layer by exchange bias 20 
after undergoing a brief annealing and in-field cooling 
process. The Nb capping layer protects the FeMn from 
oxidation and is not superconducting. 

Samples were patterned for four-terminal current-in- 
plane resistance measurements by mechanical masking 
during sputtering. The lateral dimensions of the samples 
were 4.3 mm x 1.6 mm. The critical temperatures of all 
samples were determined by ac resistance measurements 
with current of 10 fj,A, corresponding to a current density 
less than 3 x 10 5 A/m 2 , low enough to be in the linear 
response regime. T c was defined to be the temperature 
at which the resistance dropped to half its normal state 
value. 

The results for the T c measurements are summarized 
in Fig. 1. The T c of the trilayer shows a strong depen- 
dence on the superconductor thickness close to a critical 
thickness, df, where the sensitivity to ferromagnetism 
is enhanced. No superconductivity is observed above 36 
mK for d s < df « 20.5 nm. 

We have verified the magnetic configuration of our 
structures on simultaneously sputtered samples of larger 
lateral size using a SQUID magnetometer. Fig. 2 shows 
a plot of magnetization vs. applied field, H, for a sample 
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FIG. 1: Critical temperature vs. Nb thickness for a series of 
Py(8)/Nb(d s )/Py(8)/Fe 5 oMn 5 o(8)/Nb(2) samples (all thick- 
nesses are in nm). The various symbols represent different 
sputtering runs. The solid line represents the theoretical fit 
as explained in the text. Inset: R vs. T for a d a —21 .5 nm 
sample illustrating the difference between T c for the P and 
AP states. 
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FIG. 2: Magnetization vs. applied field for a d s = 23 nm 
sample measured at T = 4.2 K. At H ~ ±10 Oe the free bot- 
tom Py layer switches while the pinned top Py layer switches 
at around -500 Oe. Inset: minor loop measured at T — 4.2 K 
showing good switching of the free Py layer. 



The inset to Fig. 2 shows a minor hysteresis loop illus- 
trating that applied fields of ±100 Oe switch the trilayer 
fully between the P and AP configurations. The nearly 
zero net magnetization observed at -100 Oe suggests very 
good AP alignment, while the nearly saturated magneti- 
zation observed at ±100 Oe indicates good P alignment. 
Similarly, well-defined alignment of the P and AP states 
can be achieved at temperatures in and below the super- 
conducting transition. 

Measurements of T c p and T^ p were performed by al- 
ternating the applied field between ±100 and -100 Oe, 
while the temperature was slowly decreased through the 
transition region. The largest shift in critical tempera- 
ture, AT C = T^ p — T P , should occur in samples with the 
Nb thickness close to df. The inset to Fig. 1 shows a 
plot of R vs. T for a sample with a nominal thickness 
d s = 21.5 nm, measured in a dilution refrigerator, with 
a T c = 1.42 K. Two distinct transitions are observed for 
P and AP alignment, with a separation in temperature 
AT C « 9 mK for this case. Samples with d s ~ 22 nm 
have T c 's between 2 and 3 K and exhibit values for AT C 
of only a few mK, similar to results obtained previously 
in other F/S/F systems ,iLi2ii2i There is no observable 
difference between the P and AP state for samples with 
d s > 26 nm. 

Figure 3 shows a plot of AT C vs. T c for nine samples. 
The largest observed AT C for our Py/Nb/Py trilayers 
is about 20 mK for a sample with d s = 20.5 nm and 
T c = 0.385 K. The data are somewhat scattered for the 
thinnest Nb layers due to the increased sensitivity of T c 
to small variations of thickness and growth conditions. 
Nonetheless, our samples always show that T p < T^ p , a 
result that is opposite to what was observed by Rusanov 
et a/. 14 in similar Py/Nb/Py trilayer systems. 




with d s = 23 nm taken at 4.2 K, illustrating the typical 
spin-valve behavior of the trilayer. The narrow hystere- 
sis loop near H — shows the switching of the free Py 
layer with a coercive field H c = 5 — 10 Oe, while the 
wider loop shows switching of the pinned layer, shifted 
to nonzero H due to the exchange bias with the FeMn. 



FIG. 3: AT C vs. critical temperature for a series of 
Py(8)/Nb(d s )/Py(8)/Fe 50 Mn 5 o(8)/Nb(2) illustrating the dif- 
ference in T c between the P and AP state. The fit to the data 
is obtained using the theory of Fominov et alw— as outlined 
in the text. 
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The critical temperature of F/S/F trilayers in the P 
and AP states has been calculated theoretically by sev- 
eral groups 4i2iSiiSSiSi The usual approach involves solv- 
ing the Usadel equations in the dirty limit, which for 
the superconductor implies Is < £,bcs — ^s7/7r 2 fcsT c b , 
and for the ferromagnet lp < hvp/E exl where Is and lp 
are the electron mean free paths in S and F. Here, Tj! is 
the transition temperature of the bulk superconductor, 
vs and vp are the Fermi velocities in the S and F mate- 
rials, and 7 = 1.7811. These simplified theories do not 
consider different electronic properties for the majority 
and minority spin bands of the F material. 

We compare our data with the theory of Fominov et 
al^ The following equations, which describe the critical 
temperatures T c p and T AP for the P and AP cases, are 
obtained in the limit of a thin S layer with a constant 
superconducting gap A, and a strong ferromagnet with 
E ex > A: 



In 



T c o 
fj 



hi 



T C Q 



2 d.Tf) 
2 d s T C AP J ' 



(1) 



0, (2) 



where £g = y/hDs/2irkBT c0 and T c0 is the critical tem- 
perature for an isolated superconducting layer of thick- 
ness d s . Fominov et al. make the important point that 
the existence of a significant dependence of T c as a func- 
tion of the relative magnetization angle for d s > £s is due 
to the fact that the critical temperature of the trilayer is 
suppressed as compared to that of the isolated Nb layer, 
i.e. T c <C T c o. Consequently, the condition for which this 
theory is valid, d s <C £ = y/hDs/{2irkBT c ), is consider- 
ably weaker than the condition d s <C £s, because £ ^> £5. 
In the limit of thick ferromagnets, the tanh functions in 
refiSi are set to 1, and the functions Vh and W in Eqns. 
n and [21 become 



V h 



Ps£s 



{l-i)p F /2k h + R B A 



W = Re{V h } (3) 



where kh — y E ex / HD p is the inverse coherence length in 
F and RbA is the boundary resistance times area of the 
S/F interface, a parameter whose value reflects both the 
quality of the interface and the Fermi surface mismatch 
between the S and F materials. Eqns. (1-3) produced 
the fits to the Py/Nb/Py data shown in Fig. 1 and 3. 

Estimates for the parameters appearing in the theory 
were obtained by performing additional measurements on 
bulk and thin film samples. Since the F layer is treated in 
the thick limit and its thickness remains fixed for all our 
samples, we have used a bulk value for the resistivity of 
Py, namely pp = 123 nfi mi24 By contrast, the thickness 
of the S layer in our trilayers changes, and thus we have 
measured ps as a function of the thickness on bare Nb 
thin films. In addition, the variation of T c q with thick- 
ness was also measured on the same Nb films. The ex- 
plicit dependencies of ps and T c were taken into account 



in Eqns. ^andGI 25 The coherence length was obtained 
by performing perpendicular field measurements on the 
bare Nb films, giving £s ~ 6 nm in the thickness range 
of our data. Taking the limit of T c — * in Eq. ^ f° r 
the behavior as d s approaches the critical thickness df, 
results in the relation df/£s = 2e c |14| where C = 0.577 
is the Euler constant. Using this constraint in Eq. [31 one 
can obtain an estimate for the boundary resistance: 



R B A « 2e c p s (d c s 



e s 



(4) 



where the value of ps is taken at the critical thickness. 
After constraining RbA as shown above and using the 
measured values for the resistivities and £s> kh is the 
only remaining fit parameter. 

Using Eq. 0]with df = 20.5 nm and £g = 6 nm gives 
the estimate RbA — 1.5 ffi m 2 , which when utilized in 
Eq. yields a fit that follows the T c vs d s data very 
well, as shown in Fig. 1. The fit to the T c vs d s data is 
somewhat insensitive to the value of kh, however. That 
parameter is tightly constrained by fitting to the AT C vs 
T c data. The results of that fit are illustrated in Fig. 3, 
showing good agreement to the data with kh — 1.0 nm -1 . 
In comparison, independent estimates of kh using the 
values of vp and If discussed earlier are 0.8 and 2.0 nm -1 
for the majority and minority spin bands, respectivelysAS. 
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FIG. 4: AT C vs. critical temperature for a series of 
Ni(7)/Nb(d s )/Ni(7)/Fe 5 oMn 5 o(8)/Nb(2)p±£ The fit to the 
data is obtained using the theory of Fominov et alw— 



The excellent fit shown in Fig. 3 motivated us to apply 
the theory of Fominov et al. to our previously-reported 
data on Ni/Nb/Ni trilayersi^ The results of the fit to 
the AT C vs T c data from the Ni/Nb/Ni trilayer are il- 
lustrated in Fig. 4, and also show excellent agreement, 
even though it is not obvious a priori that the Ni lay- 
ers in those samples are in the dirty limit. The values 
pp = 33 nfhn, £ s = 6 nm and df = 16.5 nm were 
used in the fit, which gave RbA = 2.3 ffim 2 for the 
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Ni/Nb interface and kh = 0.5 nm -1 . We have made 
independent measurements of the Nb/Ni interface re- 
sistance using current-perpendicular-to-plane resistance 
measurements of Nb/Ni multilayers, and find the value 
RbA = 2.35 ±0.25 Kl m 2 , in excellent agreement with 
the value obtained from the fit to the T c vs. ds data. Our 
independent estimate of kh varies over a broad range due 
to uncertainty in determining the value of the diffusion 
constant (or mean free path) in Nii& From the measured 
resistivity, we obtain values of / p ranging between 7 and 
70 nm, depending on what value we take for the product 
PfIf for NiiSi Combining that with the values for the 
exchange energy E ex = 0.115 eV and the Fermi velocity 
vp = 0.28 x 10 6 m/sfii we obtain values for kh ranging 
from 0.16 — 0.5 nm" 1 . The value corresponding to the 
shorter lp agrees with the value from the fit to the data 
in Fig. 4. 

The question remains open as to why Rusanov et 
al*^ observe inverse spin switch behavior, T c p > T C AP , 
whereas we observe the standard behavior, T p < T^ p . 
The most obvious difference between our samples and 
theirs is that we use exchange bias to pin the magnetiza- 
tion direction of one Py layer, whereas they rely on the 
different coercivities of the two layers. But the switching 
data in their micron-scale samples show a clear plateau, 
which suggests that they have achieved a good AP mag- 
netization configuration. A second comment is that they 
observe a difference between T c p and T C AP even when the 



Nb layer is very thick, 60 nm, whereas sensitivity to the 
ferromagnet orientation is limited to our samples with 
d s < 28 nm. Variations in resistance or T c have also 
been observed in F/S bilayers due to domain formation 
during magnetization switching^Zi2& But Rusanov et al. 
state that the features indicating the inverse spin switch 
effect in their trilayers were not observed in bilayers. This 
fact, combined with their data in micron-scale samples 
that appear to be single-domain, argue against any role 
of domains in producing the inverse effect. 

In summary, we observe similar spin-switch behavior 
in Py/Nb/Py and Ni/Nb/Ni trilayers - both S/F sys- 
tems with strong ferromagnets. The results from both 
systems are fit well with the dirty-limit Usadel theory of 
Fominov et al~ This success is somewhat unexpected 
given that this oversimplified theory assumes identical 
electronic properties (density of states, Fermi velocity, 
and mean free path) for the majority and minority spin 
bands of the ferromagnetic material, whereas Py is known 
to have a strong spin-scattering asymmetry. The success 
of a dirty-limit theory in Ni is also surprising, and may 
be due partly to strong diffusive scattering of electrons 
from the S/F interfaces. 
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